The INTEGRAL source IGR J16328-4726: a High Mass X-ray 

Binary from the BeppoSAX era 



M. Fiocchi^ A. Bazzano^ A. J. Bird^ S. P. Drave^ L. Natalucci^ P. Persi \ L. Piro^-^ P. 

Ubertini ^ 

ABSTRACT 

We report on temporal and spectral analysis of the INTEGRAL fast transient 
candidate IGR J16328-4726 observed with BeppoSAX in 1998 and more recently 
with INTEGRAL. The MEGS X-ray data show a frequent micro activity typical 
of the intermediate state of Supergiant Fast X-ray Transients and a weak flare 
with duration of ~4.6 ks. The X-ray emission in the L5-10 keV energy range 
is well described through the different time intervals by an absorbed power law 
model. Comparing spectra from the lower emission level up to the peak of the 
flare, we note that while the power-law photon index was constant (~ 2), the 
absorption column density varied by a factor of up to ~6-7, reaching the value 
of ~2xl0^^ cm~^ at the peak of the flare. 

Analysis of the long-term INTEGRAL /IBIS light curve confirms and refines the 
proposed ~10.07 day period, and the derived ephemeris places the BeppoSAX 
observations away from periastron. 

Using the near and the mid-IR available observations, we constructed a spectral 
infrared distribution for the counterpart of IGR J16328-4726, allowing us to 
identify its counterpart as a High Mass OB type star, and to classify this source 
as a firm HMXB. 

Following the standard clumpy wind theory, we estimated the mass and the radius 
of the clump responsible of the flare. The obtained values of Md ~ 4 x lO^^g 
and Kci — x 10^ km are in agreement with expected values from theoretical 
predictions. 
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1. Introduction 

The X-ray fast transient source IGR J 16328-4726 was discovered by Bodaghee et al. 
(2007). Fiocchi et al. 2010, making use of the IBIS/INTEGRAL data, tentatively classified 
it as a Supergiant Fast X-ray Transient (SFXT) candidate on the basis of the transient and 
recurrent nature, the spectral properties observed during the outbursts that lasted only a few 
tens of minutes/hours, the location in the Galactic Plane, the faint quiescent emission (< 
2.5 X 10~^^ erg cm~'^) and a flux dynamic range >170 in the 0.3-10 keV (XRT/Swift) and 
>165 in the 20-50 keV (IBIS/INTEGRAL) bands. The source characteristics and behavior 
are reported in Fiocchi et al. (2010) and references therein. An additional detection was 
reported in the Swift BAT 54 month all-sky survey (Cusumano et al. 2010) with an averaged 
flux of (3.2 ± 0.8) X 10~^^ erg cm~'^ in the 15-70 keV energy band. 
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Fig. 1. — Left: Lomb-Scargle periodogram of the 18— 60keV IBIS/ISGRI light curve of 
IGR J16328-4726 showing the detection of the orbital period at 10.068±0.002 days. Right: 
IBIS/ISGRI 18— 60keV phase folded light curve using the best orbital period determination 
of 10.068 days and a zero phase ephemeris of MJD 52651.164. The locations of the BeppoSAX 
observations within orbital phase using the stated ephemeris are also indicated. 

The physical mechanism producing short X-ray flares of the SFXTs is still under dis- 
cussion, and several scenarios have been proposed. The most frequently discussed scenario 
is the clumpy stellar wind model (in't Zand 2005; Walter & Zurita Heras 2007) where dense 
clumps in an isotropic stellar wind cause short X-ray flares when they accrete onto a neutron 
star. A few SFXT systems can be explained by an eccentric NS orbit around a massive star 
with a disk-like anisotropic stellar wind (Sidoli et al. 2007). The SFXT emission could be 
explained also by a gating effect caused by the extremely strong NS magnetic field (~ 10^^ 
G) and rotational properties (Bozzo et al. 2008a). 
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Proximity to the Black Hole Candidate 4U 1630-47, has previously made any periodicity 
study based on long-term hard X-ray light curves difficult, but Corbet et al. (2010) reported 
a significant period at ~ 10 days from ICR J16328-4726. 

Here we report on analysis of 5e]?^oS'y4X observations performed in 1998, a re-analysis 
of the INTEGRAL/IBIS light curve, and an archival search for other observations that might 
allow us to confirm or discard the SFXT nature of ICR J16328-4726. 



2. INTEGRAL observations and data analysis 

The INTEGRAL dataset used included all public data available for the period up to 
30*^^ September 2010 (MJD 52650 through MJD 55469). The total on-source exposure on 
ICR J16328-4726 during this period was ~9.2Ms. Light curves for both ICR J16328-4726 
and 4U 1630-47 were extracted from science window (pointing) images constructed using 
the standard 0SA9 analysis software, in the 18-60 keV energy band. The light curves were 
then filtered to remove science windows (ScW) with an exposure of less than 200 s and/or 
in which the sources were at an off-axis angle of greater than 12 degrees as such ScWs have 
large intrinsic uncertainties in flux measurement. The flnal, flltered ICR J16328-4726 light 
curve has an on source exposure of ~7.1 Ms. We performed a timing analysis on the derived 
IGR J16328-4726 light curve using standard Lomb-Scargle techniques (Lomb 1986, Scargle 
1982) to search for periodic signals in the data, and monte-carlo randomisation to determine 
both the significance of, and uncertainty in, the derived periodicities (for further information 
on the randomisation techniques applied see Drave et al. 2010 and references thererin). 

3. INTEGRAL timing analysis 

The Lomb-Scargle periodogram of the 18-60 keV IBIS/ISGRI light curve of IGR J16328- 
4726 (Figure [H left panel) shows the detection of a period at 10.068 days, which we interpret 
as the likely orbital period. A randomisation monte-carlo determined the 1-sigma uncer- 
tainty on that period to be 10.068±0.002 days. This period is consistent, within errors, with 
that reported in Corbet et al. (2010) and provides an independent confirmation of the or- 
bital period of the binary system. The detection in that work used a SWIFT /BAT data set 
spanning MJD 53359 through MJD 55104. This data set began after the last period of high 
X-ray activity in 4U 1630-47 ceased at ~MJD 53300. The INTEGRAL/lBlS data set used 
in this work does span the earlier region of high activity of 4U 1630-47. However there is 
only a small amount of contamination from this bright source in the IGR J16328-4726 light 
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Fig. 2.— Top: the MECS light curve of IGR J16328-4726 in the 1.5-10 keV energy band 
with bin size 300 s, using BeppoSAX observations performed in 1998. Bottom: the MECS 
mosaic image in the 1.5-10 keV energy band during Phase G reported in the top panel. The 
white circle shows IGR J16328-4726 position and the 4' radius circle is the region from which 
we have extracted the spectra. The bright source is the nearby BH candidate 4U1630-47. 
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Fig. 3. — Left: the MECS spectra shown together with the power law models at the peak 
emission (Phase G) and at the low emission level (Phase L) are indicated with squares and 
stars, respectively. Middle: Contour level of the column density versus the unabsorbed flux 
in the 2-10 keV measured from the spectrum during Phase G. The contours correspond to the 
68%, 90% and 99% confidence level. Right: Contour level of the column density versus the 
unabsorbed flux in the 2-10 keV measured from the spectrum during Phase L. The contours 
correspond to the 68%, 90% and 99% confidence level. 



curve and a comparative analysis performed using only data from the same MJD range as 
Corbet et al. (2010) provided a consistent detection of the period. The higher exposure in 
the full data set allowed a more accurate determination of the period and its uncertainty 
however and thererfore it is this data set that is used in the remainder of this work. 

Folding the light curve using the best orbital period determination of 10.068 days and 
a zero phase ephemeris, taken as the first INTEGRAL observation in the light curve, of 
MJD 52651.164 yields the folded light curve shown in Figured] (right panel). This shows a 
single broad emission maximum above a quiescent level, but the statistical quality makes it 
difficult to derive further information. The shape of the profile is similar to those observed 
from other SFXTs such as ICR J17544-2619 (Clark et al. 2009) and ICR J16465-4507 (Clark 
et al. 2010). The times of the BeppoSAX observations within this orbital phase using the 
stated ephemeris are also indicated. 



4. BeppoSAX observations and data analysis 

Results presented here are from the Low- Energy Concentrator Spectrometer (LECS; 
Parmar et al. 1997) and the Medium-Energy Concentrator Spectrometer (MECS; Boella et 
al. 1997) on-board BeppoSAX. The High Pressure Gas Scintillation Proportional Counter 
(HPGSPC; Manzo et al. 1997) and the Phoswich Detection System (PDS; Frontera et al. 
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Fig. 4. — Left: The column density versus the unabsorbed flux in the 2-10 keV energy range, 
for each studied period reported in figure 1. Right: The photon index versus the unabsorbed 
flux in the 2-10 keV energy range, for each studied period reported in figure 1. The green 
points are the photon index versus the unabsorbed flux in the 20-50 keV energy range using 
the IBIS/INTEGRAL data (from Fiocchi et al. 2010). 



1997) instruments were not considered as their field of view included also emission from the 
nearby Black Hole Candidate 4U1630-47. 

The LEGS and MEGS event files and spectra were generated with the Supervised Standard 
Science Analysis (Fiore, Guainazzi & Grandi 1999). The MEGS spectra were accumulated 
in circular regions of 4' radius and publicly available response matrices were used. The off 
axis corrections have been taken into account using the apposite response matrices. 
The Point Spread Function (PSF) of the MEGS instruments has been modeled as the sum of 
2 components: a Gaussian describing mainly the detector PSF and a King Profile describing 
mainly the optics PSF (Boella et al. 1997). The final result is a fully analytic differential 
PSF with a width of 2.5' at 6.4 keV, 80% of encircled energy radius. Since the distance be- 
tween the Black Hole 4U 1630-47 and IGR J16328-4726 is ~ 14', the Black Hole contribution 
within the spectral extraction radius of the IGR J16328-4726 is widely below 1% (being this 
the systematic error added to the data following the indications given in Fiore, Guainazzi & 
Grandi (1999) for the standard spectral analysis) . 

In the archival MEGS images (1.5-10 keV) IGR J16328-4726 was detected above 5 sigma 
during three different observations performed on 1998-03-07 12:23:21.0 UTG, 1998-03-19 
14:52:11.0 UTG and 1998-03-26 17:16:48.5 UTG. The LEGS image analysis (0.5-3.5 keV) 
showed that IGR J16328-4726 was not detected above 5 sigma during these three observa- 
tions. Log of these observations are shown in Table 1, reporting start times, exposure times 
and the IGR J16328-4726 count rate for each observations. 
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5. The BeppoSAX MECS results 



Figure 2 (upper panel) shows the MECS light curve of IGR J16328-4726 with a bin size 
of 300 s. As previously observed by IBIS on board of INTEGRAL satellite (Fiocchi et al. 
2010), the source shows highly recurrent micro-activity with durations from tens of minutes 
to several hours. 

To search for spectral changes, we accumulated X-ray spectra during several time in- 
tervals, as indicated on the light curve. The X-ray spectra are all well fitted by an absorbed 
power-law in the 1.5-10 keV energy range. Details of the time intervals, and results of this 
analysis for each selected period are reported in Table 2. 
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Fig. 5. — Left panel: Color-coded 2MASS image centered on IGR J16328-4726 made from 
J(blue), H(green) and K(red) individual images. The circle is the XRT uncertainty at 99% 
confidence level from Grupe et al. 2009. Right panel: the infrared energy distribution of the 
IGR J 16328-4726 counterpart. 



Figure |3] (left panel) shows the MECS spectra and residuals with respect to the absorbed 
power law models at the peak of emission (phase G) and at the lowest emission level (phase 
L), indicated with squares and stars, respectively. A plot of the evolution of the spectral 
parameters of the source is provided in Fig. |H The power-law photon index remained fairly 
constant (see Fig. HJ right panel) including both the MECS and IBIS data from Fiocchi et 
al. 2010, whereas the absorption column density underwent changes (Fig. HJ left panel). In 
particular, we measured a large increase (factor of ~6-7) in the absorption column density 
from the lower emission level (phase L) up to the peak of the fiare (period G) when the 
column density reached a value of ~2xl0^'^ cm~^. To test the significance of the measured 
variation in the absorption column density, we report in Fig. |3] the A^h-FIux parameters 
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confidence contours for the relevant spectra in phases G and L (middle and right panel). 

The contours correspond to 68%, 90%, and 99% confidence level. A very significant increase 
in the absorption column density with the unabsorbcd fluxes is clearly detected. Although 
adding an iron line does not signiflcantly improved the flt, we considered this component to 
estimate the iron line equivalent width upper limits. With the iron line centroid flxed at 6.5 
keV, we obtained equivalent width upper limits of 440 eV and 262 eV (3(T confidence level) 
for the peak of the flare (period G) and the lower emission level (phase L), respectively. 

6. Summciry of the main outburst characteristics 

The MEGS light curve of IGR J16328-4726 (Figure 2) shows highly recurrent micro- 
activity with durations of tens of minutes and a strong outburst (its peak is reported as epoch 
G in Figure 2). During the outburst detected with the MEGS instrument, the source was 
not signiflcantly detected (> 5a) in the corresponding LEGS image, that had an exposure 
time of 775 s. 

The spectral analysis of the outburst peak was reported in Table 2 (period G) for comparison 
with other activity period of the IGR J16328-4726 and summarised in Table 3, where we 

reported start time, duration and spectral parameters. Using these values and assuming a 
distance of 3 kpc and 10 kpc, the X-ray luminosity (2-10 keV) spans from 10^^ erg s~^ to 
-|^g36 gj.g g-i^ respectively. 

To determine the best position, and as the source is highly absorbed, we extracted 
the MEGS image in the 4-10 keV energy band during its outburst (G period reported in 
Figure 2). IGR J16328-4726 was detected at (51.5 ± 3.0) x lO'^ count/s at RA(J2000) = 
16h32m38.1s and Dcc(J2000)=-47° 23' 56.8" with a error box of 35 arcsec (Figure 2, lower 
panel). Note that the circle reported in this flgure is the radius for spectral extraction and 

Table 1: BeppoSAX Observations 



start Time 


Exposure time 


Count 


s-2 




LEGS 


MECS 


LEGS 


MEGS 


UTC 


s 


s 


[0.5-3.5 keV] 


[1.5-10 keV] 


1998-03-07 12:23:21.0 


12290 


29530 


n.d." 


7.3 ± 0.2 


1998-03-19 14:52:11.0 


9912 


27786 


n.d.« 


10.0 ±0.2 


1998-03-26 17:16:48.5 


11694 


31265 


n.d ° 


2.8 ±0.1 



" n.d. means that the IGR J16328-4726 was not detected above 5 a. 
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not the position error box. 

7. Infrared identification 

The transient X-ray source IGR J16328-4726 has been identified with the near-IR source 
2MASS1632379 1-4723409 (Grupe et al. 2009), shown in Fig. (left panel). The 2MASS 
catalogue reports only the J magnitude and gives an upper limit in H and K. From the 
images we have obtained the complete photometry of this source using the calibration re- 
ported in the headers. The measured magnitudes are J=14.32±0.06, H=12.29±0.03 and 
K=11.19±0.02. The source is reported in the DENIS survey as J1632379-472341 and has 
been detected also by the Spitzer/GLIMPSE survey. The magnitudes reported in the 4 IRAC 
bands obtained from the GLIMSPE catalogue (G336. 7492+00.4223) are M3,6;„„=10.21±0.05, 
M4.5Mm=9.81±0.06, M5.8;.m=9.63±0.05 and M8.o/.m=9.52±0.04. From the observed near-IR 
colors of the source J-H=2.03 ± 0.06 and H-K=1.10 ± 0.03, we derive a visual extinction of 
Av=20. The source is not detected in the I band (DENIS survey) and it is not reported in 
other optical catalogues. Combining the near and the mid-IR observations we have obtained 
the spectral infrared distribution for IGR J16328-4726, shown in Fig. |5] (right panel) as open 
and filled circles. Using the reddening law derived by Rieke and Lebosky (1985), we have de- 
reddened the observed spectral points (triangles in Fig. 13 right panel). This spectral energy 
distribution suggests the identification of IGR J16328-4726 as a OB star. In addition the 
de-reddened NIR colors of (J-H)o=-0.11±0.06 and (H-K)o=-0.12±0.03 are compatible with 
both OB Supergiant and OB main sequence, when compared with the colors for Supergiant 
and main sequence stars reported by Ducati et al. (2001). Although IR spectroscopy of this 
source is needed to confirm this result, the X-ray behavior and infrared analysis strengthens 
the idea that this source is a SFXT candidate. 



8. Discussion 

The INTEGRAL satellite has modified our understanding of high mass binary systems, 
showing the existence of a new population of compact objects orbiting around a massive 
supergiant star and exhibiting unusual properties, being either extremely absorbed, or char- 
acterized by very short and intense fiares. 

The SFXTs are characterized by short and bright fiares, each with a duration between 10^ s 
and 10^ s, reaching io^^~^'^ erg s~^ at their peak. In these transients, both the duration and 
shape of the bright fiares are variable and the mechanism producing the transient emission 
is still not clear. X-ray spectra (0.1-100 keV) during fiares are well described either by a fiat 
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power law below 10 keV (Sidoli et al. 2006), with a photon index of 0-1 and a high energy 
cut-off around 10-30 keV, or by a bremsstrahlung model with a temperature of 15-40 keV. 
A low intensity flaring, corresponding to an intermediate state between the quiescence and 
the flare peaks, is usually present, characterised by an X-ray luminosity of IQ-^^^^^ erg s^^, 
as discovered thanks to XKT/ SWIFT monitoring campaign of SFXTs (Romano ct al. 2010, 
2011). In this state the X-ray spectrum is softer, with an absorbed power law with photon 
index of 1-2 (Sidoli et al. 2008, Romano et al. 2010). At the moment there are 11 firm 
established SFXTs and a similar number of candidates (Sidoli 2011). 

A time-rcsolvcd BeppoSAX spectral analysis of IGR J16328-4726 covering a significant 
range of orbital phases shows the source X-ray emission could be well described in the 
different selected time intervals by an absorbed power law model. While the photon index 
of the power law remained constant, the absorption column density was highly variable. We 
measured a significant rise in the A^h from ~3 to 20x10^^ cm~^ across the transition from 
the low emission level (phase L) to the peak of the flare seen by SAX(phase G). Wc note 
that this does not correspond to the assumed periastron of the system based on the peak 
of the folded IBIS light curve but is likely associated with the accretion of a well-defined 
clump of wind. Similar behavior is usually observed during the occurrence of X-ray eclipses 
in some of the sgHMXBs and SFXT (Bozzo et al. 2008, 2011). The neutron star could be 
obscured by its supergiant companion and for this reason the X-ray emission, produced by 
the accretion close to the surface of the NS, is progressively absorbed in the photosphere of 
the supergiant star. 

With this source being confirmed as an HMXB, its location is possibly associated to the 
spiral arms of our Galaxy, as usually observed with INTEGRAL HMXB population (Dean 
et al. (2005), Lutovinov et al. (2005), Bodaghee et al. (2007, 2011)). Assuming a distance 
of 3 kpc and 10 kpc for the spiral structure of the Galaxy, the X-ray luminosity (2-10 keV) 
during the flare spans from 10^^ erg s~^ to 10^^ erg s~^, respectively. During the low emission 
level these values spans from 4x10^^ erg s~^ to 4x10^^ erg s~^. The X-ray characteristics 
observed with BeppoSAX show the frequent micro activity typical of the intermediate state, 
with a luminosity value lying between the quiescence and the flare peaks, and a soft X-ray 
spectrum (F ~ 2). 

Soon after the discovery of SFXT, to explain the bright and short outbursts, in't Zand 
(2005) proposed the sudden accretion of material from the clump wind of the supergiant; 
Negueruela et al. (2006a) suggested that SFXTs orbits should be highly eccentric to explain 
the low luminosity in quiescence. Later, the idea of a clump and spherically symmetric 
wind has been associated with an eccentric and/or wide orbit (Walter & Zurita Heras 2007; 
Negueruela et al. 2007). Based on the shape of the hght curve observed during the 2007 
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Table 2: Best-fit parameters of IGR J16328-4726 during different time intervals of the Bep- 

poSAX observations (see Fig. 1) 



Interval 


OBS. Tstart 


Interval Tstart 


Tstart Phase 


Exp. 


Rates 


JVh 


r 


^unabs 


xVd.o.f. 






MJD 




(s) 


(10~^c/s) 


(10^^ cm-^) 




(10~l^erg cm-^ s~l) 




A 


1998-03-07 
12:23:21.0 


50879.563 


0.036 


1521.2 


7.5±0.2 


7.7±1.5 


2.0±0.3 


12 7+^ ° 
'^''■'-1.5 


80/105 


B 


1998-03-07 
12:23:21.0 


50879.901 


0.070 


4713.5 


7.1±0.4 


9.511:1 


2 4+0.7 


13 4+^-'' 


42/54 


C 


1998-03-07 
12:23:21.0 


50879.988 


0.079 


1582.1 


14.9±0.9 


14.9l^-« 


1 7+1-^ 
^■'-0.9 


34.6l« 


29/36 


D 


1998-03-07 
12:23:21.0 


50880.041 


0.084 


4308.9 


8.2±0.4 




9 9+0.7 


1^ !i+7.2 
10..3_3 g 


29/36 


B 


1998-03-19 
14:52:11.0 


50891.681 


0.240 


3210.9 


11.1±0.6 


lotl 


2.2±0.6 


2211 


31/53 


F 


1998-03-19 

14:52:11.0 


50891.748 


0.247 


9577.5 


7.0±0.3 




2.0±0.4 


ii-3l?:^ 


69/93 


G 


1998-03-19 
14:52:11.0 


50891.949 


0.267 


1614.9 


30.1±1.4 


18.91-73 


2.0±0.5 


1041™ 


42/67 


H 


1998-03-19 
14:62:11.0 


50891.968 


0.269 


1569.6 


14.0±0.8 




1.8±0.8 


23 4+19 


33/33 


I 


1998-03-19 
14:52:11.0 


50892.029 


0.274 


6534.1 


8.9±0.4 




2.1±0.5 


17 4+^-* 


60/85 


L 


1998-03-26 
17:16:48.5 


50898.720 


0.939 


28460.6 


2.8±0.1 


3.3±1.3 


1.7±0.3 


3.4±0.4 


70/101 



The continuum spectral model is an absorbed power law. Here, ATh is the absorption column density, T is the power-law 
photon index, Funabs is the unabsorbed flux in the 1.5-10 keV band. In the last column we report the value of the x^/d-o.f for 
each spectrum. Exp. is the good time interval for the spectral extraction. 



Table 3: Characteristics of IGR J16328-4726 during the strongest outburst 



OBS. Tstart 
UTC 


Duration 
(ks) 


Integrationa Time 
(s) 


(10=2 cm-2) 


r 


J^unabs 
(IQ-l^erg cm-2 s-1) 


ium2-iokeV 
(erg B~l) 


xVd.o.f. 


1998-03-19 14:52:11.0 


4.6 


1614.9 


18.9lti 


2.0±0.5 


104lf, 


j()35 


42/67 



The continuum spectral model is an absorbed power law. The ATh is the absorption column density, F is the power-law 
photon index, Fu^abs is the unabsorbed flux in the 1.5-10 keV band. Outburst duration is the extrapolated activity period 
assuming as beginning of the outburst the first bin during which the source was detected while the integration time is the 

good time interval for the spectral extraction. 
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outburst from the SFXT IGR Jl 1215-5952 displaying periodic outbursts, Sidoli et al. (2007) 
proposed the presence of a second denser wind component in the form of an equatorial wind 
disk from the supergiant donor, inclined with respect to the orbital plane. The flaring activity 
was explained with a wind that is inhomogeneous and anisotropic. 

All models predict that fast X-ray flares can be produced by sporadic capture and 
accretion onto the neutron stars hosted in the SFXTs. 

The wind accretion theory allow us to estimate the physical and geometrical parameters 
of the clumpy responsible of the transient emission. The duration of the main observed flare 
during the BeppoSAX observation of our source (phase G) is Atflare ~4.6 ks, a typical time 
scale on which the neutron star crosses a dense clumpy. Neglecting the orbital velocity of 
the neutron star as reporting in Lepine et al. 2008, from the duration of the flare we obtain 
the radius of the clump accreted by the neutron star: 

i?d ^ l/2t;^ Atflare (1) 

where is the relative velocity between the clump and the neutron star. 

During the time of the neutron stars periastron passage the absorption column density 
increase and it could be computed from the equation (4) given in Bozzo et al. 2011: 

Nn ^ 1.3 X lO^'vtA^^cm-' (2) 

where ciskpc is the source distance in units of 3 kpc and v^s is the wind velocity in 

10® cm s^^. 

Using the A^h measured by the spectral analysis at the peak of the flare (~1.9xl0^^ cm^^) 
and assuming a distance of 3 kpc, Eq. |2] allow us to estimate a wind velocity of ~1900 km/s. 
This value agrees with the /3-velocity law (Castor et al. 1975) which gives a typical velocity 
of 1000-2000 km/s for an OB supergiant and for a companion star losing mass in the form 
of a steady, homogeneous and spherically symmetric wind. 

Assuming this wind velocity, we found the radius of the clump accreted by the neutron star 
Rc\ ~ l/2v„Atflare — 4.4 X lO^km. In the framework of the Bondi-Hoyle-Lyttleton accretion 
theory, since only matter within a distance smaller than the accretion radius is accreted, the 
mass of the clumpy is 

Mel = (i?cl/^acc)'Macc, (3) 

where i?acc=2G'MNs/^^w — 1-0 ^ lO^km, Mns=1-4M0 and Mace is the mass accreted 
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during the flare. The X-ray flux extrapolated in 0.1-10 keV energy range during the flare 

Funabs=Lx/ {'^T^d?) = {GM^^Mg^c/ Rns) / {4:nd'^) and the continuity equation M^^c = P v^t^ -RLc 
(where d is the source distance and -Rns=10 km is the neutron star radius) allowed to estimate 
the mass accreted onto the neutron star. We derived Macc=lO~^M0 year~^ (for a distance 
of 3 kpc) and Mace ^ 2.5 x 10^^ g, thus 

Mci = (Rci/RsccfM^, ~ 6 X 10^2^ (4) 

Ducci et al. (2009) computed the expected clumpy characteristics in the framework of the 
Bondi-Hoylc accretion theory taking into account the presence of clumps. The expected 
values of M^i and i?ci according to this clumpy wind model are in agreement with the 
measured values for IGR J16328-4726. We conclude that IGR J16328-4726 can be added to 
the hst of Supergiant Fast X-ray Transients as it shows all of the characteristics typical of 
that class of object. 

The authors acknowledge the AST flnancial support via ASI/INAF grants 1/033/10/0. 
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